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We calculate the spatially resolved tunneling conductance of topological superconductors (TSCs) 
based on ferromagnetic chains, measured by means of spin-polarised scanning tunneling microscopy 
(SPSTM). Our analysis reveals novel signatures of MFs arising from the interplay of their strongly 
anisotropic spin-polarisation and the magnetisation content of the tip. We focus on the deep Yu- 
Shiba-Rusinov (YSR) limit where only YSR bound states localised in the vicinity of the adatoms 
govern the low-energy as also the topological properties of the system. Under these conditions, we 
investigate the occurence of zero/finite bias peaks (ZBPs/FBPs) for a single or two coupled TSC 
chains forming a Josephson junction. Each TSC can host up to two Majorana fermions (MFs) 
per edge if chiral symmetry is preserved. Here we retrieve the conductance for all the accessible 
configurations of the MF number of each chain. Our results illustrate innovative spin-polarisation- 
sensitive experimental routes for arresting the MFs by either restoring or splitting the ZBP in a 
predictable fashion via: i) weakly breaking chiral symmetry, e.g. by the SPSTM tip itself or by 
an external Zeeman field and ii) tuning the superconducting phase difference of the TSCs, which is 
encoded in the An- Josephson coupling of neighbouring MFs. 

PACS numbers: 74.45. +c, 74.55.+V, 74.78.-w, 85.25.-j 


I. INTRODUCTION 

Recent experiments in hybrid superconducting devices 
m have revealed strong evidence for the existence of 
Majorana fermions (MFs) [4H6], which are neutral non- 
abelian zero-energy quasiparticle excitations in many- 
particle systems. The unambiguous confirmation of their 
discovery and their successful experimental manipula¬ 
tion will be an evidence for this exotic kind of statistics 
and in addition will constitute fertile ground for develo¬ 
ping topological quantum computing [?,. The first ex¬ 
periments claiming the discovery of MFs have been per¬ 
formed in platforms based on semiconducting nanowires 
jl, with strong spin-orbit coupling (SOC), a type of ar¬ 
tificial topological superconductors (TSCs) predicted in 
Ref. |8| One of the most distinctive features of an isolated 
MF is the appearance of a zero bias peak (ZBP) of height 
2 e 2 1h in the tunneling conductance [9] 10!. However, a 
ZBP is not unique to MFs m and thus a number of 
experiments [12] have scrutinised the findings concerning 
the possible presence of MFs in such devices. 

The particular debate opened the door to new ideas for 
engineering MFs without involving semiconductors but 
by instead employing devices based on conventional SCs 
in the presence of some kind of magnetic texture. One 
of the first proposals along this direction m pinpointed 
that a magnetic chain with randomly ordered spins de¬ 
posited onto a SC can harbor MFs at its edges. Later on 
it was realised [HUS] that in the absence of SOC and 
other external fields the presence of a helical magnetic 
texture is the minimal requirement for engineering a TSC 
[6|. More importantly, TSCs from magnetic chains on su¬ 
perconducting substrates can be manipulated and probed 
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FIG. 1. Two coupled topological superconductors (TSCs) 
based on magnetic chains probed by means of spin-polarised 
scanning tunneling microscopy (SPSTM). The tip electrons 
owe a spin-polarisation along the P direction. Each chain 
can host 1 or 2 MFs per edge by virtue of chiral symmetry. 
The two TSCs are interfaced by a thin insulator yielding a 
Josephson junction since the superconducting phases on the 
two chains differ by Sip. By tuning the orientation of P, the 
value of Sip, as also by controllably violating or restoring chi¬ 
ral symmetry, one can imprint Majorana fermion (MF) sig¬ 
natures on the tunneling conductance dli/dVi. 


using spin-polarised scanning tunneling microscopy (SP¬ 
STM) [163, allowing the spatial visualisation of the MF 
wavefunctions and thus providing a more reliable method 
for diagnosing TSCs. This advantage further motivated 
a plethora of theoretical proposals involving helical mag¬ 
netism HZ]. Remarkably, it has been also shown that 
helicity (or SOC) can be artificially engineered using suit¬ 
able external fields for the case of antiferromagnetically 
ordered chains [18]. If instead SOC is present, both fer¬ 
romagnetic jT9d26j and antiferromagnetic mm chains 
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on a SC can exhibit topologically non-trivial properties 
and host MFs. As a matter of fact, the former situation 
appears to have been recently realised in the laboratory 
|2]. Nonetheless, further experiments at lower tempera¬ 
tures and with improved energy resolution are required 
for confirming the presence of MFs [2] [25]. 

In this work we illustrate a number of yet unexplored 
experimental fingerprints of MFs, part of which can be 
directly tested in the existing devices [2], thus promising 
an unambiguous identification of the MFs. Key feature 
of our approach is properly taking into account the mag¬ 
netic spin-polarisation content of the SPSTM tip, which 
could also provide an alternative explanation for the very 
low value of the measured conductance in the above men¬ 
tioned experiment. Previous studies in such systems have 
focused on a non-magnetic tip in the normal [23] [25] 
or superconducting phase m ■ Note however that spin- 
selective Andreev processes due to MFs have been previ¬ 
ously studied in Ref. [28] for nanowire-based TSCs. 

Here, following the method of Ref. jjTj we retrieve 
the tunneling conductance for ferromagnetically ordered 
chains within the microscopic Yu-Shiba-Rusinov (YSR) 
[29H32] type of model, extracted in Ref. HU The latter 
electronic bound states dominate the low-energy proper¬ 
ties of the substrate SC in the presence of the magnetic 
adatoms. If the atomic spins can be treated as classical, 
the YSR states become fully responsible for the topologi¬ 
cal properties of the hybrid device. As it has been shown 
in Ref. HU such topological YSR bands can harbor 1 or 2 
MFs per edge in the presence of chiral symmetry. 

Specifically, we investigate the tunneling conduc¬ 
tance properties and the emerging zero/finite bias peaks 
(ZBPs/FBPs) for i) isolated chains supporting either 
number of MFs per edge and ii) pairs of coupled chains 
with the same or different number of MFs per edge. Cen¬ 
tral feature of our analysis is the electronic spin-polari¬ 
sation induced by the MFs, which exhibits a strongly 
anisotropic coupling to the tip electrons, transferring 
this distinctive behaviour to the tunneling conductance. 
Evenmore, depending on the spin-polarisation of the 
magnetic tip, chiral symmetry can be weakly violated 
locally in a controlled manner. The latter can drastically 
affect the conductance spectra and the qualitative fea¬ 
tures of the ZBPs. In the case of the coupled chains we 
additionally consider a finite phase difference in the order 
parameters of the SCs, yielding 47r-periodic Josephson 
MF couplings 331 . which can further control the loca¬ 
tion of the peaks. This collection of fingerprints appear 
feasible to be looked for in existing platforms [2] and ca¬ 
pable of unveiling the presence of MFs. 


functions and the inter-MF coupling in the case of two 
tunnel-coupled TSCs are both discussed in Sec. |IIB| The 
SPSTM tip Hamiltonian and its coupling to the MFs is 
contained in Secs. II C| and |IID| Finally, in Sec. [HE 


present the steps for calculating the tunneling conduc¬ 
tance by employing the Keldysh formalism, restricted to 
the Andreev processes induced by the MFs. 


A. Topological ferromagnetic YSR chain 


In a previous work, Ref. HU three of the present au¬ 
thors extracted the microscopic topological model de¬ 
scribing the YSR bands of a ferromagnetic chain de¬ 
posited on a superconducting substrate with Rashba 
SOC. The YSR bands arise from the overlap of YSR 
bound states which lie energetically deep below the ener¬ 
gy scale set by the bulk superconducting gap. For 
extracting this low energy model, one starts from the 
Hamiltonian of the bulk SC with Rashba SOC 

Hsc = ^ E'd K kTz + ar z (k x z) ■ a - A T y a y ] 'J k , (1) 
k 


where the Pauli matrices cr and r are defined in spin 
and particle-hole space respectively, a = ( r z a x , <J y ,T z cr z ), 
while 4/]”, = ^-kp ^-kO is t ^ ie corresponding 

spinor. The operators ^ a create electrons of momen¬ 
tum k and spin-projection a. Here corresponds to 
the free electron energy dispersion, a denotes the SOC 
strength and A the bulk superconducting gap. Straight¬ 
forward manipulations allow us to obtain a Hamilto¬ 
nian describing the electrons of the SC, localised at the 
adatom sites. The corresponding Schrodinger equation 
in site space i, j has the form (see Ref. HU) 


— £(j)i (2) 

j 


with (j)\ = ( u , Vi^) and the Hamiltonian 


TTUpM 2 


( 7iv F M 2 Ty(jy - Mr z cr z ) 6 i; 


+ M 2 (Gl_ 0 T Z - GUjT.ay + Ft_ jT yCJy - F?_ jT y) 


( 3 ) 


II. THEORETICAL MODEL AND METHODS 

In the following section we present in detail our theo¬ 
retical approach and model. In Sec. EE we review the 
topological properties and model of a single ferromag¬ 
netic YSR chain. The form of the related MF wave- 


with the Fermi-level density of states (DOS) in the nor¬ 
mal phase of the SC, vp. M = JS denotes the ferromag¬ 
netic energy scale, with J the exchange energy between 
the SC electrons and the adatoms with spin S. The solu¬ 
tion of Eq. © determines the energies and wavefunctions 
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of the YSR midgap states, while the coefficients read 
G s {r ) 

7 XVp 

F s (r) 

7 XVp 

G a (r) 
iixvp 
F a (r ) 

ilXVp 


= cos (kpar/vp) sin [kp\r\ 
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= sin (kpar/vp) sin ( kp\r\ 
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4 / \ 

1 r 1 

/ 7rk F \r\ 5 

= sin (kpar/vp) cos ( kp\r\ 


/ 2 
/ nk F \r\ 

- ( 7 ) 


where £o is the coherence length of the SC, kp (vp) the 
Fermi wave-vector (velocity) and r = (i — j)a, with the 
adatom spacing a. Note that G a,s ( 0) = F a,s ( 0) = 0. The 
indices s and a denote functions which are symmetric or 
anti-symmetric under inversion r —>> — r. In the rest of 
the manuscript we use kp in units of 7r/a, a in units of 
vp , M in units of 1 /( 7 xup) and £o i n units of a. 

By transferring to k -space we obtain the Bogoliubov- 
de Gennes (BdG) Hamiltonian 


'Hk = tkT z - v k r z (jy + (A + V k )r y (jy - d k r y - Bt z (j z , ( 8 ) 

where we have introduced B = A/(7t vpJS) and 


OO 

cos (Ska) 

8=1 

with 

ts = 

2 A 

pis 

- ^8 5 

7 XVp 

(9) 

OO 

vs sin (Ska) 
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with 

V § = 

2 A 

pi a 
^8 ^ 
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( 10 ) 

OO 

Vs cos (Ska) 

8=1 

with 

F s = 

2A 

( 11 ) 

OO 

d§ sin (Ska) 

with 

d ( 5 = 

<1 ^ 

CM t? 

( 12 ) 


8=1 


source of chiral symmetry breaking is a Zeeman field ap¬ 
plied along the y axis, i.e. B y . Instead B x and B z fields 
preserve the latter. The B y field can be either applied 
globally or only near the edge, since it will always hy¬ 
bridise the 2 MFs sitting at the same edge. Notably in 
a SPSTM experiment the tip itself is magnetised, and 
depending on its magnetic orientation it can controllably 
violate or preserve n. The latter property has significant 
ramifications when measuring the tunneling conductance 
with the SPSTM technique, and can influence the pos¬ 
sible observation of the ZBP or its quantisation with a 
single or double unit of conductance ( 2 e 2 /h). 


B. Majorana wavefunctions and their coupling 

By solving Eq. ©> for a finite-size chain, one can re¬ 
trieve the MF wavefunctions which have the special form 

= > u i,l,n j u i,\,n > u i y i,n) ? ( 13 ) 

with the wavefunction components satisfying 
JT a |rq ;(J?n | 2 = 1/2. The latter normalisation im¬ 

plies the anticommutation relation for the corresponding 
MF operators {y n ,7™} = 5 n ,m- For a low energy 
description, one can focus only on the MF sector and 
neglect the rest of the BdG quasiparticles which lie above 
the bulk superconducting gap. Thus in the low energy 
limit, we can approximate the YSR state operators 
with the MF operators, i.e., fii :Cr = J2 n u i,(r, n ln- In the 
presence of suitable symmetries and ideal conditions 
(e.g. infinite chains) the MFs operators remain unpaired 
[4], while in any realistic situation one has to introduce 
also couplings for the MFs, yielding the following general 
Hamiltonian 

— 2 ^ ^ « / ^nm7n'Tm • (14) 

n,m 


The BdG Hamiltonian above resides in symmetry class 
BDI 0EQ, with time-reversal symmetry 0 = /C ( 0 2 
/). chiral symmetry U = r x and charge-conjugation 
S = ry/C, where /C denotes complex conjugation. The 
particular symmetry class supports a Z topological in¬ 
variant m allowing an integer number of MFs per chain 
edge. The detailed diagram of 0, 1 and 2 MF phases 
per edge has been extracted in Ref. [21] as also in others 
works within different frameworks [20] 35 . Topological 
YSR chains with multiple MFs edge modes can be found 
also in two-dimensional systems [36]. 

External perturbations which violate chiral symmetry 
n (or equivalently 0 ) enforce the system to reside in 
symmetry class D, which allows up to a single MF per 
edge. This implies that if the TSC resides in the topo¬ 
logical phase with 2 MFs per edge, the application of an 
infinitesimally weak n-violating perturbation m, will un¬ 
avoidably hybridise the 2 MFs, splitting them into finite 
energy (proportional to m) bound states. However, the 
1 MF per edge phase remains unaffected by such a weak 
field. Symmetry analysis demonstrates that the simplest 


Solely residing on the above MF description is a good 
approximation only as long as the coupling elements 
Mnm are much smaller compared to the YSR bandstruc- 
ture gap. As long as this is case, the matrix elements can 
arise due to the following reasons: i) a weak chiral sym¬ 
metry breaking term, e.g. m, coupling for instance 2 
MFs at the same edge previously protected by n, ii) fi¬ 
nite size of the chains which can allow the overlap of MFs 
primarily located at the far edges, giving rise to quasi¬ 
particles with finite energy splitting Se and iii) coupling 
M of neighbouring MFs located at the edge of two dif¬ 
ferent tunnel-coupled TSC chains. For the purposes of 
our discussion MF coupling matrix elements of the first 
type will be discussed at a phenomenological level, while 
for the others we will explicitly calculate the couplings 
by employing the BdG formalism and considering a par¬ 
ticular model for inter-chain tunneling, respectively. 

For inferring the coupling of edge MFs due to the inter¬ 
chain tunneling, we focus on the electronic degrees of the 
two substrate SCs and consider that they are separated 
by a thin insulating film yielding a Josephson junction 
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with a corresponding superconducting phase difference 
Sip. The latter can be imposed by inducing a supercur¬ 
rent flow through the junction or by gluing together the 
very left and right edges of a single chain in order to form 
a ring through which we can thread flux. To this end, we 
assume that the electrons of the two substrate SCs are 
coupled via the Hamiltonian 


^ = E + h.c. 

iJ 


(15) 


In particular, here we consider the profile 


T — t 

-•-I,] I* 


1 - sgn (i-j) 


exp 


| ai — bj\ — (a + b ) 

l 


(16) 


In the above we have already assumed that one TSC 
chain resides on the positive sites i, j > 0 and the other 
on the negative ones. In this manner, the projector [1 — 
sgn(i • j)]/2 allows only interchain tunneling, while we 
have also assumed that the tunneling strength decays 
exponentially with the distance \ia — jb\ over a characte¬ 
ristic decay length /, with a, b denoting the adatom spa- 
cings of each chain. Moreover, in the above Hamiltonian 
we have incorporated the spatially varying superconduc¬ 
ting phase profile pi, which here is assumed to have the 
form cpi = sgn(i)Sip/2. The corresponding MF matrix 
elements that one obtains take the following form 

M nm = 41m E ■ (17) 


C. SPSTM tip Hamiltonian 


where p^ — p^ocP and v a = 1. We now introduce the 

spin polarisation degree P s = ranging from — 1 

to 1. The latter values occur for fully spin-polarised tips. 
P s can strongly vary depending on the tip material, while 
it is in principle possible to achieve complete polarisation 
(P s = ± 1 ) by fabricating the tip using a half-metal, for 
which one of the spin bands does not cross the Fermi 
level (for more details see Ref. 37). 

For the rest of our analysis, it is eligible to perform 
a rotation and diagonalise the tip Hamiltonian in spin 
space using the property 

P • <y = Pa z R , R = exp {irja y /2) ex.p(i'da z /2) . ( 21 ) 
Thus the tip Hamiltonian becomes 

77tip, i = E " e ^)4,a,A, CT ,i • ( 22 ) 

k,cr=db 

with a = dz labelling the two eigenstates of a z with eigen¬ 
values ±1 and e*, ?cr = e;, — aP. 


D. Coupling between the SPSTM tip and the MFs 

The SPSTM tip originally couples to the electronic 
density of the superconducting substrate at site i, via 
the tunneling Hamiltonian 

'H'riP-MF.i = E ( T m4,ct,;Vv + H.c.) . (23) 

k,cr 

After diagonalising the tip Hamiltonian in spin space we 
obtain 


For the purposes of this work we model the SPSTM tip 
at site i, as a lead of spinful electrons under the influence 
of a spin splitting field P. Moreover, here we assume 
that the tip feels a voltage V*, which is responsible for 
driving the coupled system out of equilibrium and leads 
to the tunneling current. Previous works on the tun¬ 
neling conductance for such systems have considered a 
non-magnetic tip in the normal [23] [25] or the supercon¬ 
ducting phase m- The tip Hamiltonian has the form 

"HtIP.j = E ~ eV i) 5 a0 ~ P ' ^]4,a,t c lk,^,i ’( 18 ) 
k,a,(3 

Here by keeping the index i for the tip electrons and the 
voltage we allow our formalism to address the case of 
a multi-tip SPSTM, while in the most general case we 
should include an index i also to the spin-polarisation P. 
The latter can be parametrised in the following way 

P = P(cos # sin t?, sin d sin 77 , cos rf ). (19) 

The presence of P modifies the DOS for the spin up and 
down electrons, p a {P)i at energy E. In particular the 
normalised DOS read 

V a = Pa(E)/p(E) with p{E) = E P<r(E), (20) 


■HtIP-MF,* = E ( V k,a,i,n 5 k,a,i ~ 7 n , (24) 

k,cr,n 

with the matrix elements 

Vk,cr,i,n — ^ ^ T k,iRq,a ,r ^ j i,a',n • (25) 


E. Tunneling conductance 


After having set the stage for calculating the tunneling 
conductance, we can proceed with deriving the Heisen¬ 
berg operator for the current through the tip located over 
site i 


Ii(t) = -eN T ip(t) = “ e E ^ (' 

k,a 




= “T- E ^TIP-MF,*, C' k a i C k 
k,a 


it). 


(26) 


The current operator in the Schrodinger picture reads 


h 


E + K,a,i,n%a,i) 7 n ■ (27) 

k,cr,n 
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The expectation value of the current operator Ii(t) = 
(. Ii(t )) is calculated from the expression 

h(t) = iE Im [ V k,a,i,n (7 n(t)c k , a ,i(t))] 


k,cr,n 

2e 


= -ir E Re VUi.nGk.w&t) > (28) 


k,a,n 


which involves the lesser mixed Green’s function 


G k,a,i,n( f y *') = «(7n(i')cfe,«T,i(i)) • (29) 

For calculating the current one can employ the Keldysh 
formalism and introduce the respective Keldysh-contour- 
ordered Green’s functions. After following this route we 
find that for retrieving the final expression for the current 
we will need the retarded MF Green’s functions 


C(*. t') = ~ lQ (t - t') . (30) 

Since here we are interested in the non-equilibrium steady 
state, one can follow the method of Ref.flOland show that 
the expectation value of the current operator is given by 
the formula 

+00 

h = j dco Ti(co) [n F (u} - eVi) - n F (-co + eV))] , (31) 

— OO 


with the Fermi-Dirac distribution np(uj) at energy uj and 
the transmission coefficient 


Ti(u) = Tr 


g K ( w )r(- u )g A Hf*( w ) 


(32) 


according to terminology of the Landauer-Biittiker for¬ 
malism. Note however that in contrast to the usual Lan- 
dauer formula for ballistic trasport in normal conductors, 
the present formula involves one electron and one hole 
Fermi-Dirac distribution. This is a direct consequence 
of the involvement of MFs which are equal superposi¬ 
tions of electrons and holes. The matrix elements for 
the linewidth hermitian matrices T 1 (uj) are given by the 
expression 

= 27I~ ^ ~ e k,cr) • (33) 

fc,cr 

As in Ref. uni we adopt the wideband approximation ac¬ 
cording to which the linewidth matrix elements can be 
considered to be energy independent. In particular, we 
assume that T*,^ = T (the k independence reflects the 
wideband approximation) and set the DOS of the tip 
to be approximately p^Ep)- The latter consideration 
yields 

KM = Km = H * 1 + T* Ui ’ m ’ (34) 

with r = 2np{E F )\T\ 2 , u[ n = and the 

unit vector P = P/P. Note that for convenience we set 
T = 1, a convention which we will follow throughout the 


remainder of the manuscript. For retrieving the trans¬ 
mission coefficent one has to calculate the retarded and 
advanced MF matrix Green’s functions given by 

§ R (co) = (wi-iM + I Re r) , (35) 

where Q a (uj) = [Q r (uj)]^. In the above we have conside¬ 
red that the self-energies of the MF matrix Green’s func¬ 
tions contain only the linewidth functions, in the spirit 
of the wideband approximation as in Ref. 10. 

For the rest of our discussion, we will focus on the 
zero temperature tunneling conductance given by the fol¬ 
lowing expression 

c\T 2e 2 

wr^ %(eVh (36> 


III. RESULTS AND DISCUSSION 


In the following paragraphs we present our results con¬ 
cerning the tunneling conductance for a single or two cou¬ 
pled TSC magnetic chains. In particular, the case of a 
single chain is examined in Secs. GHZ] an d |IIIB| where we 
analyse TSCs with 1 and 2 MFs per edge, respectively. 
Later on, we consider the situation of two coupled TSC 
chains with the superconducting phases of the substrate 
SCs differing by Sep, thus allowing a 47r-periodic Joseph- 
son coupling between the neighbouring edge MFs. Specif¬ 
ically, in Sec. m. C]we consider two chains, each of which 
harbors a single MF per edge. In Sec. IIIDpIIE ) we in¬ 
vestigate Josephson junctions consisting of a TSC chain 
with a single MF per edge and a TSC chain with 2 MFs 
per edge, where we assume that the SPSTM tip couples 
to the chain hosting 1(2) MFs per edge. In Sec. 


IIIF 


considered a Josephson junction of two TSC chains each 
of which hosts 2 MFs per edge. In the whole discussion 
particular emphasis will be given to chiral symmetry vio¬ 
lation and restoration, and its implications for observing 
the ZBPs which are considered smoking gun signatures 
of MFs. 


A. One TSC chain with 1 MF per edge 


In the present paragraph we discuss the case of an iso¬ 
lated chain with a single MF per edge as shown in Fig.[2j 
Similarly to the result of Ref. [10] the tunneling conduc¬ 
tance in this case takes the simple Lorentzian form 

dh _ 2e 2 (Itj 2 

dV h ( eV) 2 + (ri a )2’ 1 ^ 

from which one recovers the ZBP with one quantum of 
conductance. In the above and from now on, we set 
Vi = V for convenience. The broadening is given by 
the expression 


r l = r?y T 

L aa L u 'z a 


1 + PsP ' O’ 


Ui 


(38) 
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FIG. 2. TSC based on a ferromagnetically ordered atomic 
chain probed by SPSTM. The tip is considered to couple lo¬ 
cally, i.e. only to the electronic density directly below it. For 
long chains only the y a MF is seen by the probe, while for 
short ones both MFs contribute. 


Crucial feature of our analysis is the inclusion of the 
magnetic characteristics of the SPSTM tip which opens 
new perspectives for detecting MFs. In fact, one ob¬ 
serves that if P a u^ a = — u^ a the linewidth term be¬ 
comes T l aa = V]Tu\ a Ui, a - If the tip becomes fully spin 
polarised, so that = 0, the tunneling conductance will 
also go to zero and the ZBP will disappear from the tun¬ 
neling spectra. Essentially, when these conditions are 
met, the spin polarisation of the tip-electrons is antipa¬ 
rallel to the electronic edge polarisation ( a) i = u\ a cr u^ a 
induced by the MF. Therefore the tip-MF coupling be¬ 
comes zero since tunneling between the tip and the sub¬ 
strate electrons is spin conserving and cannot take place 
for antiparallel polarisations. 

As shown in Fig. [3] the MFs of both sides of a sin¬ 
gle chain induce an electronic spin-polarisation which is 
confined in the xz plane, as a result of the complex conju¬ 
gation symmetry 0 = JC which forces the spin-part of the 
wavefunction to be real. Note that a similar spin polari¬ 
sation profile was previously retrieved for nanowire-based 
TSCs in Ref. 38j as a result of the common features of the 
two models. The particular distinctive feature, allows us 
to employ a SPSTM tip for unveiling the MFs. In partic¬ 
ular, one obtains a characteristic anisotropic dependence 
of the tunneling conductance on the angles ($, rf) which 
determine the orientation of the tip-magnetisation, as it 
has been also pointed out in Ref. |28] within a different 
context. 

In order to make a connection to the realistic situa¬ 
tion, potentially relevant to the experiment of Ref. [2j we 
further take into account the influence of the remaining 
MF away from the tip. If the overlap of the MF wave- 
functions is non-negligible due to the short length of the 
chain, a finite coupling of the form i5e 7 a 7 e will appear 
leading to finite energy excitations. More importantly, 
for short chains the tip generally couples to both MFs. 
In this case we have 

( _° c o ) and r ‘ = ( (rfj* r\ e e ) > (39) 



Site i 



FIG. 3. SPSTM tip coupled to a single MF for the case of 
a sufficiently long chain of TV = 100 sites, (a) Spatial pro¬ 
file of the tunneling conductance for a spin unpolarised tip 
(P s = 0). (b) Tunneling conductance at site i = 1 for a tip 
polarised along the x axis for varying P s . The modification 
of the polarisation degree alters the profile broadening, (c) 
Local spectral weight (\ui\ 2 ) and electronic spin polarisation 
((<x)) of the two MFs. Note that the spin polarisation along 
the y axis is exactly zero as a result of © = JC symmetry, (d, 
e) Tunneling conductance at site i — 1 as a function of the 
angles ($, rf) for the case of a fully polarised tip P s = d=l. The 
present maps reveal the strong electronic spin anisotropy in¬ 
duced by the MF y a and constitutes one of its characteristic 
signatures. Parameters: £o — 80, kF = 6.0, a m 0.01 and 
M = 0.85. 


yielding the modified tunneling conductance formula 


dh 2e 2 


12 [det(Re P) + Se 2 


+(eC) 2 (It,) +(r‘ e ) +2(Re 2 ri e -Im 2 ri e 


detP 


(40) 


)]} 


{[' 


n 2 


(eV) 2 - de 2 - det(Re P) + (eV) 2 (P aa + P ee ) 


M = 


-1 
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FIG. 4. (a) Spectral weight of MF wavefunctions and their in¬ 
duced electronic spin-polarisation in the case of a short chain 
consisting of N — 15 sites. The two MFs hybridise and form 
finite energy quasiparticle excitations, (b) Matrix elements 
of the linewidth matrix Y l . For each case corresponding to a 
fully spin-polarised tip, we present only the non-zero elements. 
Parameters: £o — 80, Jgf — 6.0, a — 0.01 and M — 0.85. 


For V = 0 one obtains the simple result 


dh 

dV 


= 2 


2e 2 


r 2 - ir* | 

aa ee I x ae I 


V=0 


h ri a q e - Re 2 ri e 


+ Se 2 ’ 


(41) 


which provides the ZBP height in this general case. Stri¬ 
kingly, when both MFs are accessed by the tip, the ZBP 
persists , though with a conductance with reduced spec¬ 
tral weight from the ideal value. This ZBP appears due 
the coupling of the tip to the y e MF, and directly disap¬ 
pears if we set Y\ e = Y l ae = 0. Apart from the residual 
spectral weight for IF = 0, the conductance shows finite 
bias peaks 


el/ = iy^e 2 + det(Ref i ), (42) 

as shown in Fig. [5] The tunneling conductance for these 
voltages reads 


dh 

dV 


FBPs 


2e 2 (ri a +P ee ) 2 -4Im 2 q 

h (XL + Ue ) 2 


(43) 


Notably the FBPs have a height equal to 2 e 2 /h only if 
Imf^g = 0. Otherwise, the height is lower. As inferred 




Site i 



Site i 



Site i 


FIG. 5. Spatial tunneling conductance profile in the case of 
a short chain consisting of A = 15 sites, (a) Profile at site 
i — 1 for spin unpolarised (o) and fully polarised tip along 
the x,y,z axes. For a spin unpolarised tip one finds resid¬ 
ual spectral weight for V = 0. Due to the spin anisotropy 
in the MF induced spin-polarisation we observe a significant 
difference in the height of the FBPs. When the polarisation is 
confined in the xy plane or the tip is unpolarised, the conduc¬ 
tance is practically equal to 2e 2 //i, while for a tip magnetised 
along the z axis one obtain a reduced height. We present the 
complete spatial profile for a fully spin-polarised tip along the 
x,y,z axis in (b-d). Note the spectral weight asymmetry in 
b). Parameters: £o — 80, kF — 6.0, a — 0.01 and M — 0.85. 


by Fig. s we obtain an almost quantised conductance 
for all the cases except when probing with a magnetic 
tip with polarisation along the z axis. In this case, the 
tunneling conductance is much weaker. As a matter of 
fact, this is exactly the configuration employed in the ex¬ 
periment of Ref. [2] and if our YSR model is applicable, 
our findings can provide one possible explanation to the 
highly reduced signal, aside from the unavoidable tem¬ 
perature broadening. According to our theory, orienting 
the spin-polarisation of the tip along the x, y axes will 
drastically increase the conductance value to almost a 
single quantum, while we additionally obtain that the 
ZBP can be found even for very short chains of N = 15 
sites. These remarkable findings are depicted in more de¬ 
tail in Fig. ([5| and demonstrate once again that SPSTM 
is a powerful tool for detecting these MF fingerprints at¬ 
tributed to the anisotropic MF spin-polarisation. 

Finally, if the coupling of the tip to y e is completely 
negligible, then the tunneling conductance formula reads 
(see also Ref. uni) 

dh _ 2e 2 (el/) 2 (rj Q ) 2 

dV h [(eV) 2 — Se 2 ] 2 + (eV) 2 (r^) 2 ' 


From the above expression we see that the ZBP disap¬ 
pears and splits into two FBPs appearing for eV = =b Se. 





































































B. One TSC chain with 2 MFs per edge 


In this paragraph we examine the case of a single TSC 
magnetic chain where due to the preservation of chiral 
symmetry 2 MFs appear per edge as in Fig. [ 6 ] Although 
phases with two MFs per edge have not been experi¬ 
mentally demonstrated yet, they appear prominent to be 
realised in the near future. In fact, one can engineer a 
2MF per edge phase starting from a IMF per edge phase, 
similar to the one discovered in Ref. [ 2 j As prescribed in 
Ref. HU the latter can be achieved via a topological phase 
transition effected by varying a set of paramaters such as 
the adatom spacing, the SOC strength or the applied 
magnetic field. 

Here we focus on a TSC magnetic chain in the 2 MF 
phase with the two edges being infinitely separated, al¬ 
lowing us to restrict ourselves to the MF subspace of y a 
and 7 ^. Essentially we exclude coupling between MFs of 
different edges and at the same time assume negligible 
coupling between the tip electrons and the MFs located 
at the right edge. Consequently, both MFs on the left 
edge couple to the tip, while they can also couple to each 
other with a matrix element m, arising due to weakly 
broken chiral symmetry. As we have already underlined 
the latter violation can be a consequence of the SPSTM 
tip itself, if the polarisation contains a component along 
the y axis. Under these conditions we have 


M = 


0 m 
—m 0 


and Y l 


(K d Y r; 


■ ad 
i 

dd 


, (45) 


which are identical to the ones found in Eq. (39) with 
the correspondence m Se and d e. Note that here 
the off-diagonal elements for the linewidth functions are 
crucial and cannot be neglected a priori, since the two 
chiral symmetry protected MFs have spectral weight at 
the same region. This is in stark contrast to Sec. Ella] 


also previous studies m where the overlap of the 2 MFs 
involved can be completely neglected in the infinite sepa- 



FIG. 6. TSC harboring 2 MFs per edge protected by chiral 
symmetry. For sufficiently long chains solely the two MFs 
below the SPSTM probe tip need to be taken into account 
for inferring the tunneling spectra. When owing a polarisation 
component along the y axis, the magnetic tip itself can break 
locally chiral symmetry and hybridise the two MFs, even in 
the absence of other symmetry breaking fields. 


ration limit. The tunneling conductance can be obtained 
from Eq. (41) after performing the replacement Se —>• m 
and e d. For V = 0 one obtains 


dU 

dV 


= 2 • 

v=o 


2e 2 

h r L r dd - R e 2 ri ad + m2 


(46) 


As previously, the ZBP still persists and here its height 
can be controlled in general by the chiral symmetry brea¬ 
king field m as also ImT^. As shown in Fig. [7^-b), the 
electronic spin polarisation of the MF wavefunctions is 
also in this case confined to the xz plane. In addition, 
both wavefunctions are real or imaginary. Therefore, by 
employing a tip with polarisation along the y axis one 
simultaneously induces finite values for m and ImT^. 
For a spin-polarisation of the tip in the xz plane or an 
unpolarised tip, m = ImT^ = 0. Therefore, in the latter 
cases we obtain a ZBP with double unit of conductance 
[TO] , as if the 2 MFs were unpaired [4]- 
When chiral symmetry is preserved (m = Im Y l ad = 0), 
one obtains the following profile for a general voltage bias 


dh 

dV 


2e^ 

h 


(47) 


2 det 2 (f' i ) + (eV) 2 (r * 0 ) 2 + (P ee ) 2 + 2{T\ e f 


n 2 


(eVT-detn +(eV) 2 (ri a + ri e y 


One observes that only the first summand is respon¬ 
sible for the double unit of conductance ZBP, while the 
second term of the above equation contributes beyond a 
crossover voltage where the sharp spike profile switches 
to a broad hump feature. Essentially, for very small vol¬ 
tages the two MFs behave as unpaired and beyond the 
crossover voltage they couple giving rise to two FBPs. 
This is obvious from the contribution of the second sum¬ 
mand depicted with the orange line in Fig. [7 Jd) . 

Note finally that there also special cases in which the 
presence of 2 MFs per edge can be even masked and 
misinterpreted as a single MF per edge. In fact, if chi¬ 
ral symmetry is preserved and for some particular values 
the condition r* ^ = (rjj 2 = Re 2 r ad additionally 
holds, then we obtain the expression for the tunneling 
conductance 

dh 2e 2 (r i a + iyy 

dV h (eR)2 + (ri a + rjj 2 ’ 1 ; 

which is identical to the one for a single MF per edge 
but with an effective broadening Y l aa + Y dd . The special 
condition satisfied above implies that essentially only one 
MF of the chiral symmetry protected MF pair is seen by 
the SPSTM tip. 


C. Two TSC chains with 1 MF per edge 

We now consider two coupled chains, each of which 
can harbor a single MF per edge, a situation depicted in 
Fig-! The particular setup can be useful for indirectly 
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Site i 
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Site i 


FIG. 7. (a) MF spectral weight \ui\ 2 and electronic spin- 

polarisation. The spin-polarisation has only x , z components, 
(b) Profile of the double unit of conductance ZBP at i = 1 
when chiral symmetry is preserved. There exists a crossover 
voltage at which the spike-like profile switches to a broad 
hump-feature. This reflects the contribution of two different 
sources to the conductance. The blue curve shows the full 
conductance, whereas the orange curve shows only the con¬ 
tribution of the second summand in Eq. ( [48] ) . In the inset we 
zoom around V = 0. (c) Spatial tunneling conductance profile 
for broken chiral symmetry (m = 0.1) and spin unpolarised 
tip. Parameters: N — 100, £o = 80, kF = 6.0, a = 0.01 and 
M = 0.85. 


probing the 47r-periodic Josephson effect. By assuming 
that the two chains are sufficiently long so that the MFs 
away from the junction can be excluded from our analy- 


y 










Ml 2 


-\-Sp/2 


FIG. 8. Josephson junction of two tunnel-coupled TSC chains 
each of which harbors 1 MF per edge. For sufficiently long 
chains only the two MFs near the junction need to be taken 
into account. Here only y a couples to the SPSTM tip. The 
coupling M of 7 a and 7& has form of 47r-periodic Josephson 
term. 



FIG. 9. (a) Spatial profile of the tunneling conductance for 
Sp — 0. One obtain two FBPs at eV = ±M(5p) (b) 2tt- 
periodic dependence of the tunneling conductance (i = 1) on 
the superconducting phase difference Sp. For 5p = tt one 
recovers the ZBP. We used a tunneling constant t — 0.1 and 
a decay length l = 10 in units of a = b = 1. 


sis, we obtain 


M = 


0 M \ 

-M 0 J 


and f l 




(49) 


In this case the MF coupling matrix element originates 
from interchain tunneling and has the form 


M = UmJ2 u t 




± i,3 c ' 




(50) 


Similarly to Eq. (44) one obtains for the particular case 


dh _ 2e 2 (eC) 2 (7 Q ) 2 

dV h [{eV) 2 — M 2 ] 2 + (eV) 2 (r * a ) 2 ' 


Therefore, by imposing a difference Sep between the 
phases of the two SCs, we can tune M and modify the 
location of the FBPs. The tunneling conductance is a 
27r-periodic function of Sip, resulting from the 47r-periodic 
Josephson coupling oc cos(5(p/2 ) between the y a and 7 &. 


D. Two coupled TSC chains: one with 1 MF per 
edge (below the tip) and one with 2 MFs per edge 

As previously we concentrate on the MFs near the 
junction. In this case only y a couples to the tip and 
we have 



( 0 

Mab 

M ac 1 

1 ^ / r «« 

0 

0 

M = 

-M ab 

0 

TTlbc 

,r*= 0 

0 

0 


V -M ac 

-m bc 

0 J 

1 \ 0 

0 

0 


The matrix elements M^ and M ac originate from in¬ 
terchain tunneling while m^ c originates from chiral sym¬ 
metry breaking only in the left or even in both chains. 
Weak violation of chiral symmetry does not introduce 
any modification to the wavefunction of the y a MF. Thus 
in the particular case, we obtain the conductance formula 

dh = M _ [(e]/) 2 -my 2 ( r U 2 _ 

dV h {eV) 2 [{eV) 2 _ M 2J2 + (rL) 2 [{eV) 2 _ m 2j2 - 

( 53 ) 
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FIG. 10. Josephson junction of two tunnel-coupled TSC 
chains harboring a different MF number per edge (two and 
one respectively). For sufficiently long chains, only the three 
MFs near the junction need to be taken into account. Here 
only 7 a couples to the SPSTM tip. The couplings M of y a 
to 7 b and 7 C have the form of 47 r-periodic Josephson terms. 
Chiral symmetry breaking can further mix j b and 7 C . 

with M = y/Kb + MZa + 7 tt| c . For V = 0 we observe 
that we obtain the ZBP. Indeed, this confirms the rule 
derived in Ref. |10] according to which the tunneling con¬ 
ductance of an odd number of coupled MFs demonstrate 
the ZBP. Moreover there are also two FBPs of 2e 2 /h at 
eV = However, if chiral symmetry is preserved, i.e. 

rribc = 0, we obtain 

dh = 2 e 2 (eC ) 2 ( 7 0 ) 2 

dV h [( e v)2 - M 2 } 2 + (eV) 2 (T i aa ) 2 ' 

Essentially if chiral symmetry is preserved the system 
behaves as only 2 MFs become coupled, with an effec¬ 
tive coupling M = \J M 2 h + M^ c . This becomes transpa¬ 
rent by writing ij a (M ab ^ b + M ac 7 c ) = iM^ a (M ab ~f b /M + 
Mac'Jc / M ). The orthogonal linear combination of MFs 
M ac ^ b — M ab ~f c remains unpaired [3] and unseen by the 
SPSTM tip. Therefore switching on and off the chiral 
symmetry breaking field can controllably make the ZBP 
appear or disappear providing a smoking gun signature 
of MFs in these chains. On the other hand, one by con¬ 
trolling Sip can shift the two split peaks at eV = =bM, 
adding another experimental knob for detecting MFs. 


FIG. 11 . Josephson junction of two tunnel-coupled TSC 
chains harboring a different MF number per edge (one and 
two respectively). For sufficiently long chains, only the three 
MFs near the junction need to be taken into account. Here 
only the chiral symmetry protected MF pair y a and 7 d couple 
to the SPSTM tip. The couplings M of 7 *, to 7 a and 7 d have 
the form of 47 r-periodic Josephson terms. If chiral symmetry 
is broken y a and 7 d can also have a finite coupling. 


The matrix elements M ab and M db originate from inter¬ 
chain tunneling while m ad originates from chiral sym¬ 
metry breaking in the particular or even both chains. 
Weakly breaking chiral symmetry does not affect the 
wavefunction of the single edge MF. 

Here the resulting tunneling conductance expression 
is lengthy and therefore we present it in the Appendix. 
However, one can directly infer that the ZBP persists both 
in the presence or absence of chiral symmetry. The reason 
is that the MFs below the tip become always indirectly 
coupled via their interaction with the tip, and thus the 
presence or not of chiral symmetry is only important for 
ensuring the presence of the 2 MFs. Since all three MFs 
are coupled and their number is odd the appearance of 
the ZBP was expected according to the rule of Ref. flOj 
If the couplings to 7 b become zero, then we return to the 
case of 2 MFs at the edge of a single chain discussed in 
Sec. |IIIB[ in which case the existence of a ZBP crucially 
depends on the persistence of chiral symmetry. 


F. Two TSC chains: Both with 2 MFs per edge 


E. Two coupled TSC chains: one with 1 MF per 
edge and one with 2 MFs per edge (below the tip) 


In the present paragraph we consider the second pos¬ 
sible configuration for two chains with uneven edge MF 
number, as shown in the Fig. IF In the particular case 
one obtains the matrices 


/ 0 ^ad M ab \ ^ 

M = -m ad 0 M db , r 

\ -Mab ~M db 0 J 


Ka Kd 0 \ 

(bJ* T dd 0 . 

0 00 / 

( 55 ) 


For completing our analysis we proceed with examining 
the case in which both chains harbor 2 MFs per edge and 
are tunnel-coupled. By focusing only on the MFs near 
the junction depicted in Fig. [l] we can write 


M = 


(r 



0 

m ad 

M ac 

M a b \ 



■m ad 

0 

M dc 

M db 



M ac 

-M dc 

0 

m cb 

1 

V- 

M ab 

-M db 

-m cb 

0 J 


■n 

aa 

pi 

1 ad 

0 0\ 




i \* 
add 

1 dd 

0 0 j 

M = 

( M ac 

M ab 

0 

0 

0 0 

, 1V1 

\M dc 

Mdb 

0 

0 

0 0/ 





(56) 


.(57) 
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The expression that one obtains for the tunneling con¬ 
ductance is lengthy and therefore we do not include it in 
this manuscript. However it is easy to retain the expres¬ 
sion for the tunneling conductance at h = 0, which is 
sufficient for providing information about the qualitative 
characteristics of the system. The zero voltage conduc¬ 
tance reads 


2e 2 2 m 2 cb det T l 

h f \ ^ ^ 

(det M - m a drn cb J + m 2 cb det(ReT*) 

(58) 

It is straightforward to observe that we obtain a ZBP 
but the conductance does not take a quantised value. At 
this point we can further focus on special cases. If chiral 
symmetry is at least restored for the chain which does 
not couple to the SPSTM tip (i.e. m c b = 0), the ZBP va¬ 
nishes, implying that effectively an even number of MFs 
out of the four are essentially probed by the tip. On the 
other hand, if chiral symmetry is only restored for the 
chain probed by the tip (i.e. m a d = 0 and P y = 0), the 
ZBP remains but with modified and still non-quantised 
conductance value. Note however that if m a d = 0 and 
det M = 0 then we obtain a peak of a double unit of 
conductance as if only two unpaired MFs appear in the 
system. Finally, if only detM = 0, the tunnel coupling 
matrix has a zero eigenvalue implying that the tip does 
not see the second chain and we return to the case of 2 
MFs probed simultaneously by the tip. 


dh 

dV 


IV. TUNNELING CONDUCTANCE BEYOND 
THE MF-INDUCED ANDREEV PROCESSES 


Throughout the whole manuscript the tunneling con¬ 
ductance was computed based on the approximation that 
the YSR operators can be replaced with the MF opera¬ 
tors, i.e., W = J2 n u i,o r,n7n- Essentially, we kept only 
MF-induced Andreev processes. In fact, these are the 
only possible Andreev processes which can occur, since 
the SPSTM is considered non-superconducting in the 
present discussion (see e.g. 39). The remaining contribu¬ 
tion to the conductance can only arise by single electron 
processes. Due to the bulk gap in the energy spectrum 
of the TSC chain, single electron tunneling is suppressed 
and can only occur if some YSR states become unoc¬ 
cupied, so that the tip-electrons can tunnel into them. 
This can happen for finite temperatures or/and due to 
inelastic scatering with collective modes in the substrate 
such as phonons. For more details see Ref. 1271 where 
it has been also shown that the inelastic processes are 
thermally activated. Therefore the applicability of our 
results is restricted only to very low temperatures where 
the Andreev approximation is well justified, while in ad¬ 
dition it is required for the SPSTM tip to be located very 
close to the substrate. 


V. CONCLUSIONS AND PERSPECTIVES 


We explored new distinctive MF features which can 
be measured via spin-polarised scanning tunneling mi¬ 
croscopy in hybrid devices consisting of ferromagnetic 
chains on top of spin-orbit coupled superconductors. For 
the calculations we adopted a microscopic model descri¬ 
bing Yu-Shiba-Rusinov chains which can harbor 1 or 2 
MFs per edge if chiral symmetry is present. 

For an isolated topological YSR chain with a single 
MF per edge, we showed that the tunneling conductance 
delicately depends on the direction along which the tip is 
spin-polarised. If fact, for a fully spin-polarised tip there 
can be special angles of the tip polarisation for which the 
conductance vanishes. In addition, in the case of short 
chains where MFs on both edges contribute, one finds 
that depending on the tip polarisation the signal can be 
extremely weak while at the same time other directions 
support an almost quantised ZBP. This can be relevant 
for the recent measurements of Ref. [2] which were charac¬ 
terised by a weak signal. According to our findings, the 
polarisation and Zeeman field of the experiment (i.e. z 
axis) were oriented along the direction for which the MF 
signal appears to be the lowest. Within our scenario, 
altering the field angle could lead to signal enhancement 
and allow the possible observation of the long-sought-for 
ZBP. 

Moreover, we showed that in the case of a single chain 
with 2 symmetry protected MFs per edge, the tip (or an 
additional Zeeman field along the y axis) can induce a 
weak chiral symmetry violating term which can control- 
lably modify the tunneling spectra. In addition, for the 
case of tunnel-coupled topological chains, one can induce 
a difference between the phases of the two superconduc¬ 
tors in order to modify the location of emerging finite bias 
peaks via the 47r-Josephson coupling. In fact, the tun¬ 
neling conductance could be used as an indirect probe 
of the latter. Furthermore, for two chains with different 
edge MF number, tunable chiral symmetry violation and 
restoration can switch on and off the ZBP which is a 
robust MF feature. 

Note that in experiments based on self-assembled mag¬ 
netic chains, as in Ref. ,2 . junctions can be currently dif¬ 
ficult to fabricate. Nonetheless, a number of the above 
mentioned Josephson effects can be still accessed by ei¬ 
ther inducing a supercurrent flow along a short TSC mag¬ 
netic chain (see e.g. Refs. [151 and l40l) or employing a chain 
in a ring geometry through which flux can be threaded. 
In both configurations the MFs of the left and right edges 
are coupled and feel different superconducting phases, 
thus behaving similarly to the neighbouring MFs of two 
tunnel-coupled chains. 

Conclusively, this novel MF characteristics relying on 
the MF spin-polarisation, which were extracted from a 
realistic YSR microscopic model for these chains, reveal 
new experimental methods for unambiguously detecting 
MFs in the near future. 
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Appendix: Expression for the tunneling conductance for the case presented in Sec. |III E| 


dli _ 2e 2 N{eV) 
dV ~ ~h D(eV) ’ 

where we have introduced the denominator 

DM = w 6 + W 4 [(rjj 2 + (r dd ) 2 + 2 Re 2 r ad - 2 M 2 } 

+ u 2 {[r^r*. , d - Re 2 r ad + m 2 } 2 - 2(it a + r dd ) [r aa M 2 db + r dd M 2 ab - 2(Re r ad )M ab M db \} 
+ l^aa^db + PMC _ 2(ReT l ad )M ab M db } 2 , 
and the nominator 


(A.l) 


(A.2) 


N M = {CUT + (U) 2 + 2 [Re 2 r ad - im 2 r ad ]} + 2 W 2 [r aa r dd - Re 2 r a 
-2u 2 { [(r*J 2 + Re 2 r ad - Im 2 Vj M 2 b + [(M) 2 + Re 2 - Im 2 UJ M 2 , 


+ m 2 ad ] [r aa r dd - |rMl 2 ] 

2(K a + r dd ){Rer ad )M a bM db } 


+ [^aa^db + F dd M 2 b ~ 2(ReT l ad )M a bM d b} 2 , 


(A.3) 


with M = Vn^J+M^+M 2 b . 








